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Abstract 

The effect of a random transverse field (RTF) on the wetting and layering tran- 
sitions of a spin-1/2 Ising model, in the presence of bulk and surface fields, is studied 
within an effective field theory by using the differential operator technique. Indeed, the 
dependencies of the wetting temperature and wetting transverse field on the probability 
of the presence of a transverse field are established. For specific values of the surface 
field we show the existence of a critical probability p c above which wetting and layering 
transitions disappear. 
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1 Introduction 



Recently much attention has been directed to study the wetting and layering transitions of 
magnetic surfaces Ising model. Experimental studies have motivated much theoretical works 
in order to understand and explain the growth of thin layers from only single atoms. A simple 
lattice-gas model with layering transitions has been introduced and studied in the mean field 
approximation by de Oliveira and Griffiths [1]. Multilayer films adsorbed on attractive sub- 
strates may exhibit a variety of possible phase transitions, as has been reviewed by Pandit et 
al. [2], Nightingale et al. [3], Patrykiejew et al. [4] and Ebner et al. [5-8]. One type of transi- 
tions is the layering transitions, in which the thickness of a solid film increases discontinuously 
by one layer as the pressure is increased. Such transitions have been observed in a variety of 
systems including for example 4 He [9,10] and ethylene [11,12] adsorbed on graphite. Ebner and 
Saam [13] carried out Monte Carlo simulations of such a lattice gas model. Huse [14] applied 
renormalization group technique to this model. It allowed the study of the effects on an atomic 
scale in the adsorbed layers. The lattice gas models applied to the wetting phenomena was 
reviewed by Dietrich [15]. The effect of finite size on such transitions has been studied, in a 
thin film confined between parallel planes or walls, by Nakanishi and Fisher [16] using mean 
field theory. 

The model of transverse field was originally introduced by de Gennes [17] for hydrogen-bonded 
ferroelectrics such as KH2PO4. Since then, this model has been applied to several physical sys- 
tems, like DyV0 2 , and studied by a variety of sophisticated techniques [18-21]. The technique 
of differential operator introduced by Kaneyoshi [22] is as simple as the mean field method 
and uses a generalised but approximate Callen relation derived by Sa Barreto and Fittipaldi 
[23]. The system has a finite transition temperature, which can be decreased by increasing the 
transverse field to a critical value £l c . The effect of a transverse field on the critical behaviour 
and the magnetisation curves was studied [18-21] and by Kaneyoshi et al. [24,25]. Using the 
perturbative theory, Harris et al. [26] have studied the layering transitions at T = in the pres- 
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ence of a transverse field. Benyoussef and Ez-Zahraouy have studied the layering transitions of 
Ising model thin films using a real space renormalization group [27], and transfer matrix [28] 
methods. 

On the other hand, the random systems have been known to be dominated by rare regions. 
This effect is particularly pronounced for random quantum systems at low or zero temperature 
far from critical points. Indeed, Griffiths [29] showed that the free energy is a non analytic 
function, because of rare regions. Having found all the derivatives being finite, Harris [30] con- 
cluded that this effect was very weak for classical systems. The simplest of all random quantum 
systems is the random transverse Ising model [31,32] (and references therein). 
Using the mean field theory, we have studied in a previous work [33], the wetting and layering 
transitions of a spin— 1/2 Ising model in the presence of a uniform transverse field. Our aim in 
this work is to study the effect of a random transverse field (RTF), on the wetting and layering 
transitions of a spin— 1/2 Ising system, within an effective field theory (EFT) by using the 
differential operator technique. The outline of this paper is as follows. In Section 2 we present 
the formalism and the method. In Section 3 we investigate and discuss the phase diagrams. 

2 Model and method 

Transverse and longitudinal magnetic fields are applied to a system with N coupled ferromag- 
netic layers. The Hamiltonian can be written as 

n = - £ J*s;ss - + HiS?) (i) 

<i,j> i 

where the first summation is carried out only over nearest-neighbour pairs of spins, S'f , (// = 
x,z) are the Pauli matrices of a spin— 1/2 and = J is the exchange interaction assumed to 
be constant. Hi is the longitudinal field applied on the site i, assumed to be uniform in a layer 
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k, and denned by: 

H + H 8l for k = 1 

// for 1 < k < X (2) 

H + H s2 for k = N 

where the surface fields H sl and H s2 = —H sl are applied on the first layer k — 1 and the last 
layer k = N, respectively. H is the longitudinal field. The transverse field Q is governed by the 
probability distribution law: 

ViSl) =p5{Q-Q ) + (1 -p)S{n) (3) 

with < p < 1. The case p = corresponds to the absence of the transverse field, whereas 
p = 1 is a situation with a uniform transverse field Q = Qq. 

We will use an effective field theory which neglects correlation between spins [35-36], but it takes 
into account the relations such as < (s^) 2 >= 1 exactly, where < ... > denotes the thermal 
average. This method has been used by several authors to study quantum systems [18,22-25], 
and disordered systems [18,19,21,37,38]. It is based on a cluster comprising a single selected 
site labelled % and the neighbouring sites with which it directly interacts. Hence, the part of 
the Hamiltonian containing the site i, is given by 

Hi = (£j ij s; + H i )s? + nsf (4) 

the summation runs over nearest neighbour sites j of the site i. The diagonalization of the 



operator Hi leads to the eigen values = ±^/x 2 + f2 2 , with %i = Hj^i JijSj + H t and J^- = J. 
In the next we will use the notation of = 2Sf. Following Sa Barreto and Fittipaldi [23] we 
write the approximate relation: 

„ „ , rr( g fexp(-l«,)) 

<K " >c>=< Tr(eM-m,)) " <5> 
where (5 = l/(ksT), < ... > c indicates the mean value of of for a given configuration c of all 
other spins, < ... > denotes the average over all spin configurations. 



Equations (5) are not exact for an Ising system with a transverse field, nevertheless, they have 
been accepted as a reasonable starting point in many studies [22,23]. After averaging over the 
probability distribution p, the longitudinal magnetisation of a layer k, is given by 

m z k =p< f n (x k ) > +(1 -p)< f (x k ) > (6) 

with the functions fn(x) and fo(x) are defined, respectively, by 

fn ( x ) = -^==tanhpV^T^ (7) 

f (x) =tanh(3x (8) 
Introducing the differential operator defined by: 

exp (aD)f^x) = f^x + a), = Q, 0) (9) 

where D = ^, we can write the longitudinal magnetisation of a plane k as 

m z k =p[exp(< x k > D)]fn(x)\ x=0 } + (1 -p)[exp(< x k > D)}f (x)\ x=0 \ (10) 



For a cubic lattice model, each site is in interaction with 6 neighbours, one can write the 
magnetisation ml as 

m% = p[(cosh JD+ < a z k > sinh JX>) 4 (cosh JD+ < a z k _ x > sinh JD)(cosh JD+ < a z k+l > sinh JD) 
{exp{H k D))f n {x)\ x=0 } + (1 — p)[(cosh JD+ <a z k > sinh JD) 4 (cosh JD+ < o%_ x > sinh JD) 
(cosh JD+ < a k+1 > sinh JD)(exp(H k D)) f ( 

(11) 

with the boundary conditions m z N+1 = m$ = 0. 
Neglecting the correlations, so that: 

< OiO*j...Oi >=< °t >< o] > ... < of > (12) 
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and using the differential operator, the final expression of the magnetisation ml is: 

mi = B (p, O) + B 1 (p, n)<a z k > +B 2 (p, Q) < a% > 2 +B 3 (p, Q) < a% > 3 +B A (p, Q) < a% > 4 
+ < al_ x > (B 5 (p, Q) + B 6 (p, Q) < a% > +B 7 (p, fl) < a' k > 2 +B 8 (p, Q) < a z k > 3 
+B 9 (p, Q) < a% > 4 )+ < a% +1 > (B w (p, Q) + B n (p, Q)<a z k > +B l2 (p, O) < o% > 2 
+B 13 (p, Q) < o% > 3 +B u (p, O) < o% > 4 )+ < a%_ x >< a z k+1 > (B l5 (p, Q) + B 16 (p, Q) < a% > 
+B l7 (p, n) < a z k > 2 +B l8 (p, O) < a z k > 3 +B 19 (p, Q) < a z k > 4 ) 

(13) 

where B t (p, SI) = pAi(Q) + (1 -p)A t (0), 1 = 0,1, 19, and the coefficients Ai(fi), (// = Q, 0), I = 
0, 1, 19; are given in the Appendix A. 

Similarly, the transverse magnetisation < a k >, of a layer k, can be formulated using a function 
g^x) instead of the function f^x), but in this work, we are essentially interested in the wetting 
phenomena which is essentially governed by the longitudinal magnetisation < a k >. 

3 Results and discussion 

In this section we present phase diagrams and longitudinal magnetisations of the model (1) 
described in section 2. The ground state of this model was established in a previous work 
[33]. In the following we limit our calculations to N = 20 layers. The results are found to be 
similar for an arbitrary larger number of layers iV > 20. We will distinguish between two cases, 
namely: 

i) p — 1.0: uniform transverse field, 

ii) < p < 1: random transverse field. 

The start point is a situation where all the spins are down. For fixed values of H sl /J and Q/J, 
there exist a temperature T w /J such that: 

i) for T / J < T w /J the spins of all layers are down for H/J < and up once H/J > 0, with 
the coexistence of the two cases at H/J = 0, 

ii) for T / J > T w /J and increasing the bulk field H/J the spins of the first layer will flip and 
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become up: this is the surface transition. Increasing the bulk field more and more, the spins 
of the second layer will flip up, and so on. For a sufficiently large number of layers N, the 
complete wetting is reached when the number of layering transitions is close to the number of 
layers at H/J = 0. T w /J is the wetting temperature. Similarly, the wetting transverse field 
Q w / J is obtained for fixed values of H s i/ J and T/J, and increasing the bulk field H/J. 
Hereafter, the surface field effect on the parameters T w / J and Vt w / J will be discussed for p — 1. 
The notation \ k O N ~ k is a situation where k layers are spin-up while N-k layers are spin-down 
from the top surface k — 1 to the surface of the bottom k = N. In order to study the behaviour 
of the system under the effect of the temperature, for fixed surface field H sl /J and transverse 
field, we plot in Fig. 1 the corresponding phase diagrams: (a) in absence of transverse field 
Vl/ J = 0.0 and (b) in the presence of a transverse field Vt/ J = 1.0. One can note that the effect 
of increasing the transverse field is to decrease T w /J. When increasing the temperature the 
wetting transverse field Q w /J decreases. This is shown in Fig. 2(a) for a lower temperature 
T/J = 0.5, and Fig. 2(b) for higher temperature T/J = 2.0. We can note the absence, in Figs. 
1(b) and 2(b), of the transitions: N <-> 1^. Indeed, these transitions subsist for Q/J = 0.0 
and lower temperature, Figs. 1(a) and 2(a) respectively, in agreement with the ground state 
[33]. 

Examining the magnetisation behaviour we found that it decreases rapidly in the transition 
zone either when increasing the temperature as it is shown in Fig. 3(a) for Q/J= 1.0, or 
when increasing the transverse field as it is plotted in Fig. 3(b) for T/J = 0.02. Positive 
magnetisations, in these figures, are calculated for a bulk field H/ J — > + , whereas the negative 
values correspond to H/J — > 0~. 

In order to study the effect of a random transverse field (0 < p < 1) on the wetting and 
layering transitions, we introduce the probability law given in Eq.(3). The corresponding phase 
diagrams in the space (T/J, Q/J,p) , are plotted in Fig. 4, for H sl / J = 0.95. It is important 
to note here the existence of a critical value pc, of the probability p, above which the wetting 
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transitions disappears. The numerical value corresponding to Fig. 4 isp c = 0.0885. This critical 
probability p c depends on the surface field (H s i/ J). As it is illustrated in Fig. 5, p c decreases 
linearly for values of H sl /J < 0.85, and almost linearly for H sl /J > 0.85, with a discontinuity 
point for H sl / J = 0.85. A similar result was found by Harris [34], for a disordered system with 
bond and site dilutions, showing that the critical transverse field at zero temperature presents 
a discontinuity as the concentration passes through the critical percolation concentration. 
Below, we will consider the effect of the surface field, and the probability p, on the wetting 
temperature T w /J and the wetting transverse field fl w /J. Indeed, T w /J vanishes as well as 
increasing the surface field at fixed probability value, Fig. 6(a), or increasing the probability at 
fixed surface field value, Fig. 6(b). While fl w / J does not exhibit the same behaviour. However, 
it decreases when increasing the surface field for a fixed probability value p, Fig. 7(a). Whereas 
it diverges for a fixed surface field value for sufficiently small value of p as it is illustrated in Fig. 
7(b). Furthermore, these figures show that the surface field value making T w /J = increases 
for decreasing p, whereas the surface field value for which Q w /J = is close to 1. 

4 Conclusion 

Within the effective field theory (EFT) and using the differential operator technique, we have 
studied the phase diagrams of wetting and layering transitions of a spin— 1/2 Ising model, for a 
random transverse field (RTF). It is found that, in the pure case p — 1, this system exhibits the 
same behaviour as the mean field study [33]. The dependency of the wetting temperature and 
wetting transverse field on the surface field and the probability of the presence of a transverse 
field, was investigated. We have showed the existence of a critical probability p c above which 
the wetting and layering transitions disappear. 
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Figure Captions 



Figure 1: Phase diagram in the (H/ J,T/ J) plane for a pure system (p=1.0) with the surface 
fields (H sl /J = 0.998, H s2 /J = -H al /J) a) Vl/J = 0.0, b) Q/J = 1.0. 

Figure 2: Phase diagram in the (H/J,fl/J) plane for a pure system (p=1.0) with (H sl / J = 

0.998, H s2 /J = -H sl /J) a) T/J = 0.5, b) T/J = 2.0. 

Figure 3: Magnetisations for {H sl / J = 0.95, H s2 / J = — H sl / J): 

a) as function of the temperature for VI/ J — 1.0, 

b) as function of the transverse field for T/J = 0.02. 

Figure 4: Phase diagrams in the (T/J,Q/J,p) space for (H s i/J = 0.95, H s2 / J = —H s i/J). 
Figure 5: The critical probability p c dependence on the surface field H sl /J for T/J = 0.02. 
Figure 6: Wetting temperature T w / J as a function of: 

a) surface field H s i/J for Q/ J — 1.0 and several values of the probability p 

b) probability p for Q/J = 2.0 and several values of the surface field H sl /J. 
Figure 7: Wetting transverse field Vt w / J as a function of : 

a) surface field H s \/J for T/J = 0.02 and several values of the probability p 

b) probability p for T/J = 2.0 and several values of the surface field H sl / J. 
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A Expression of the coefficients A^fi) 

Using the functions f^x), (/i = Q, 0) and the field H iy defined in the body text, the coefficients 
A k (fj,)(fx = fl, 0), k — 0, 19 for a layer i, are as follows : 



AM = (1/2) 6 x (f,(QJ + Hi) + 6^(4J + Hi) + 15/„(2 J + H t ) + 
20f^Hi) + 15/„(-2 J + Hi) + 6/„(-4 J + if,) + / M (-6J + if,)) 

^(/i) = 4 x (1/2) 6 x (/„(6J + #*) + 4/ M (4J + if,) + 5^(2J + #,)- 
5/„(-2 J + Hi) - 4/„(-4J + if,) - U-6J + if,)) 

A 2 (//) = 6 x (1/2) 6 x (/ M (6J + if,) + 2/^(4J + Hi) - 2/ M (2J + H,)- 
*U(Hi) - /„(-2 J + i^) + 2/ M (-4 J + if,) + /„(-6 J + #,)) 

As(Ai) = 4 x (1/2) 6 x (/„(6J + if,) - 3/ M (2J + if,) + 3/„(-2 J + #,)- 
f^(-6J + Hi)) 

AM = (1/2) 6 x (/ M (6J + #,) - 2/ M (4J + Hi) - fn(2J + 
4/ M (^) - / M (-2 J + - 2/ M (-4J + Hi) + fn(—6J + if*)) 

A 5 (/i) = (1/2) 6 x (/ M (6J + if,) + 4/ M (4J + ifO + 5/ M (2J + Hi)- 
5/ M (-2 J + ify - 4/„(-4 J + if,) - / M (-6J + #,)) 

A 6 (/x) = 4 x (1/2) 6 x (/„(6J + #<) + 2/ M (4J + H t ) - f,{2J + ify- 
4/^) - 5/„(-2 J + + 2/„(-4 J + if,) + fn(—6J + if*)) 

A 7 (/i) = 6 x (1/2) 6 x (/ M (6J + ifO - 3/^(2 J + Hi) + 3/„(-2 J + #,)- 
U(-6J + Hi)) 

AM = 4 x (1/2) 6 x (/ M (6J + #,) - 2/ M (4J + H t ) - / M (2 J + 
4/ M (^) - U-2J + if*) - 2fn(—4J + if,) + fii(—6J + if,)) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 
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AM = (1/2) 6 x (/ M (6J + Hi) - 4/ M (4J + Hi) + 5/ M (2 J + 
5/„(-2 J + ify + 4/ M (-4J + ify - / M (-6J + Hi)) 

AM = (1/2) 6 x (/„(6 J + if,) + 4/ M (4J + ify + 5/ M (2J + 
5/„(-2 J + Hi) - 4/„(-4 J + - / M (-6J + 

4iW = 4 x (1/2) 6 x (/„(6 J + Hi) + 2/ M (4J + - /„(2 J + 
*U(Hi) - M-2J + if,) + 2/ jU (— 4 J + if,) + /„(-6 J + ify) 

A 12 (/i) = 6 x (1/2) 6 x (/ M (6J + Hi) - 3f,(2J + H t ) + 3/„(-2 J + H t )- 
U-6J + H,)) 

AM = 4 x (1/2) 6 x (/ M (6J + - 2/ M (4J + H t ) - / M (2 J + 
4/ M (^i) - M-2J + Hi) - 2/ M (-4 J + ify + / M (-6J + Hi)) 

AM = (1/2) 6 x (/„(6J + if*) - 4/ M (4J + Hi) + 5/ M (2J + H t )- 
5/„(-2 J + + 4/„(-4 J + if,) - / M (-6 J + if*)) 

A 15 (/i) = (1/2) 6 x (/ M (6J + if*) + 2/ M (4J + Hi) - f,{2J + Hi)- 
Af^(Hi) - /„(-2 J + H t ) + 2/„(-4J + if*) + / M (-6J + 

A 16 (/i) = 4 x (1/2) 6 x (/ M (6J + if,) - 3/ M (2J + if,) + 3/„(-2 J + H t )- 
f^-6J + Hi)) 

AM = 6 x (1/2) 6 x (/ M (6 J + Hi) - 2/ M (4J + - f,(2 J + + 
^U(Hi) - U-2J + Hi) - 2/ /i (— 4J + ^) + /„(-6 J + 

A 18 (/i) = 4 x (1/2) 6 x (/ M (6J + ify - 4/ M (4J + Hi) + 5/ M (2J + H t )- 
5/ M (-2 J + if,) + 4/„(-4 J + if,) - /„(-6J + 

A 19 (/x) = (1/2) 6 x (/„(6 J + Hi) - 6/ M (4J + + 15/^(27 + H t )- 
20 f„(Hi) + 15/„(-2 J + if,) - 6/ M (-4J + tf.) + / M (-6J + if,)) 
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